Linear oscillations of a compressible hemispherical bubble 

on a solid substrate 
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The linear natural and forced oscillations of a compressible hemispherical bubble on a solid sub- 
strate are under theoretical consideration. The contact line dynamics is taken into account by 
application of the Hocking condition, which eventually leads to nontrivial interaction of the shape 
and volume oscillations. Resonant phenomena, mostly pronounced for the bubble with the fixed 
contact line or with the fixed contact angle, are found out. A double resonance, where independent 
of the Hocking parameter an unbounded growth of the amplitude occurs, is detected. The limit- 
ing case of weakly compressible bubble is studied. The general criterion identifying whether the 
compressibility of a bubble can be neglected is obtained. 

PACS numbers: 47.55.dd, 47.55.dr, 46.40.-f 



I. INTRODUCTION 

The dynamics of bubbles and droplets is of great in- 
terest for their numerous applications. Bubbly fluids are 
widely used as displacing fluids in petroleum industry. 
Small bubbles and droplets allow to efflciently control 
heat and mass transfer in heat-exchangers and reactors 
and to intensify mixing in microdevices»ii2i Whereas oscil- 
lations of drops suspended in a fluid ambient away from 
the boundaries have been scrutinized for over a century,'^ 
oscillations of drops and bubbles in contact with solid 
surfaces have only received attention for the last few 
decades. Understanding fundamental aspects of drops 
and bubbles interaction with the solid surface is closely 
related to the problem of wetting.'*'^ This knowledge is 
of practical importance because many technological pro- 
cesses deal with spreading of a liquid (a paint, a lubricant, 
or a dye) over solid surfaces. From the theoretical point of 
view, the presence of a solid surface often meets another 
problem, the contact line dynamics, which is currently 
far from being fully understood. 

Qualitatively, oscillations of a liquid drop of averaged 
radius R can be characterized by three time scales: the 
viscous relaxation time Ty = R^/v, the capillary time 
scale Tc \J pR^ jo, which is related to the period of 
the shape oscillations^ and the time scale of the acoustic 
oscillations Tq = Rjc. Here, a is the surface tension, 
p, V, and c are the the density, kinematic viscosity, and 
speed of sound, respectively. Quite often, these three 
time scales relate to each other as 



Ta < Tc < Ty. 



(1) 



This hierarchy of the times allows for the description of 
the drop oscillations within the model of inviscid incom- 



pressible fluid, where the dynamics of the drop corre- 
sponds to the shape oscillations. For instance, these in- 
equalities hold for water drops of i? > 10^^ cm. For a 
liquid drop immersed in another liquid the hierarchy of 
the time scales is the same. The only difference is in the 
meaning of p, v, and c, which should now be considered 
as some effective quantities, e.g., the mean values for the 
two liquids. 

For a gaseous bubble in an ambient liquid the situation 
is similar to the case of the drop, but an additional, the 
so-called "breathing mode," appearsii This mode corre- 
sponds to the radial (volume) oscillations of the bubble 
and is caused by the bubble compressibility. Note that 
since the densities of the gas and liquid differ consider- 
ably, the gas compressibility is non-negligible even at sub- 
acoustic frequencies, iO\,Ta <C 1, where uj-^ is the frequency 
of the breathing mode. In this case, the gas pressure in 
the bubble very quickly adjusts to the instant volume 
of the bubble. Although the pressure field in the bubble 
changes in time, its instant distribution can be considered 
as spatially homogeneous. If dissipation is insignificant 
during a period of oscillation, the gas can be d escribed by 
the adiabatic law. For this situation- Wf, = yJi^yP^JpR? , 
where 7 is the adiabatic exponent and Pg is the equi- 
librium pressure in the bubble (for a correction of uj, 
caused by surface tension see Sec.|lTl. Detailed analyses 
of damping effects can be found in Refs. [l|H. Particularly, 
these studies figure out the inequalities 



iObTy > 1, LUbTt > 1, UJbTa < 1, 



(2) 
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which can be applied to neglect the effects caused by 
viscosity, heat transfer, and acoustic irradiation, respec- 
tively. Here, we introduce the characteristic time of ther- 
mal diffusion Tf = R^ /xrn with Xm being the maximal 
from heat diffusivities of the fluid x and the gas Xg- 
it follows from ^ , the predominant damping mechanism 
for an air bubble in water is heat dissipation. This mech- 
anism can be neglected provided that R > 10~^ cm. 

If the frequencies of the shape and volume oscillations 
become comparable, the oscillations of these types start 



2 



to interact. For instance, parametric excitation of the 
shape oscihations on top of the forced breathing mode 
has been addressed by Mei and Zhou?^ who have found 
an instabihty of the radial oscillations and performed a 
weakly nonlinear analysis. Since then, the problem of the 
parametric instability has received much attention, see, 
e.g., a recent review by Feng and Leal.™ Another ex- 
ample of the interaction has been provided by Longuet- 
Higgins.— He focuses on bubble oscillations in a liquid 
and shows that the nonlinear coupling of the shape os- 
cillations can lead to the excitation of the volume mode. 
Nonlinear oscillations of an incompressible drop with the 
accurate account for the dynamics of the ambient gas 
have been studied in Ref. llJ. The nonlinear coupling of 
the shape oscillations results in generation of sound. It is 
worth noting that in the previous studies the interaction 
of the shape and volume oscillations arises as a nonlinear 
effect. In our paper we report on another, pure linear, 
mechanism of coupling of the oscillations of the different 
kinds. This mechanism of linear coupling is caused by 
the contact line dynamics. 

Last years have witnessed growing interest in under- 
standing the contact line dynamics. Although the steady 
motion of the contact line has been well studied,~ii^ there 
has been no rigorous theory for unsteady motion yet. 
What is typically applied in this situation is a simplified 
approach. The thin viscous boundary layer is neglected 
and a phcnomenological boundary condition imposed on 
the apparent contact angle is applied instead. Such con- 
dition has been proposed by Hocking for small oscilla- 
tions of the contact linet^ The velocity of the contact 
line is assumed to be proportional to the deviation of the 
contact angle from its equilibrium value (for simplicity, 
the equilibrium contact angle is considered to be 7r/2): 

|=.An.VC. (3) 

Here, C is the deviation of the interface from its equi- 
librium position, n is the external normal to the solid 
surface. The coefhcient A has the dimension of velocity 
and is referred to as the wetting or the Hocking parame- 
ter. The mostly studied cases correspond to the simplest 
limiting situations of either the fixed contact line (C = 0, 
the pinned-end edge condition) or the fixed contact an- 
gle (n • VC = 0, the free-end edge condition). Except 
for these two particular cases, the Hocking condition ([3]) 
leads to energy dissipation at the contact line. A simple 
generalization of Eq. ^ that is based on experimental 
observations'* and accounts fo r hy steresis of the contact 
line has been provided in Ref. llSl 

A consideration of natural oscillations of a hemispher- 
ical drop for the fixed contact angle has shown*^ that 
the eigenfunctions coincide with the even modes of the 
natural oscillations of a spherical drop. The problem of 
axisymmetric oscillations of a drop with the fixed con- 
tact line has been addressed both experimentally^^'*® and 
numericallyj *^i^°i^*i^^ Experimental studie o^'^i*® focus on 
the forced oscillations. In particular, eigenfrequencies are 



determined for different equilibrium contact angles. Nu- 
merical investigations deal with natural oscillations of 
an inviscid dropfi^ and with naturaP° and forced^iiS^ 
nonlinear oscillations of a viscous drop. Longitudinal vi- 
brations have recently been studied theoretically^"^ and 
experimentallyi^ In Ref. the horizontal and vertical 
orientation of the substrate is considered. For the vertical 
orientation, gravity leads to asymmetry of oscillations. 
The study by Noblin et al^ has manifested the non- 
trivial dynamics of the contact line: at a relatively small 
amplitude of vibration the contact line remains pinned, 
while at higher amplitudes it starts to move in the stick- 
slip regime. Asymmetric vibrations as well as possible 
microfiuidic applications have been discussed in Ref. H^. 

Recently, the focus of attention has been on the impact 
of the contact line dynamics on the natural and forced 
oscillations of an incompressible hemispherical drop on a 
solid substratci^i^i Axisymmetrical modes of the natural 
oscillations caused by transversal vibrations of the sub- 
strate are studied in Ref. [H. Another paper— addresses 
the nonaxisymmetrical modes of the natural oscillations 
and the forced oscillations for the longitudinal vibrations 
of the substrate, where inertia of the ambient fluid is 
taken into account. To some extent, these analyses can 
be applied to describe oscillations of an incompressible 
bubble in a liquid. However, the bubble compressibility, 
a principal feature of the present paper, can become of 
crucial importance, which has been beyond the scope of 
the previous research i^^'^^ A recent experimental study^ 
has indicated an interesting crossover, where beyond a 
certain threshold of vibration acceleration the bubble can 
split into smaller parts. Despite noticeable progress, the 
dynamics of a compressible bubble on a vibrated sub- 
strate has not been fully understood. 

In the present paper we address the behavior of a com- 
pressible hemispherical bubble on a solid substrate. The 
paper is outlined as follows. In Sec. [H] we analyze nat- 
ural oscillations. Section |llT] deals with the forced oscil- 
lations for the normally vibrated substrate. A transition 
to the case of a weakly compressible bubble is performed 
in Sec. IIVI This is the situation, when the frequency of 
the volume oscillations is high compared with that for 
the shape oscillations. Particularly, we obtain a general 
criterion identifying whether compressibility of a bubble 
can be neglected. In Sec. |V] we discuss the results and 
summarize the most important conclusions. 



II. NATURAL OSCILLATIONS 

Consider natural oscillations of a compressible gaseous 
bubble that sits on a solid substrate and is surrounded 
by a liquid. Fig. [TJ Suppose that in equilibrium the bub- 
ble is hemispherical, i.e., the equilibrium contact angle 
equals tt/2 (although this assumption is of no crucial im- 
portance, it considerably simplifies the forthcoming anal- 
ysis) . We assume that the gas density is small enough so 
that the frequency of the volume oscillations is compara- 
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FIG. 1: The geometry of the natural oscillations problem. 



ble with that of the shape oscillations. Next we impose 
the frequency restrictions ^ , which ensure insignificance 
of the damping effects caused by acoustic irradiation, vis- 
cous and heat dissipation. Finally, we admit that the 
bubble is sufficiently small so that the hydrostatic dif- 
ference of the pressure is negligible compared with the 
pressure contribution caused by surface tension. This as- 
sumption allows us to neglect gravity and ensures that 
the averaged bubble surface is hemispherical to very high 
accuracy. Formally, this situation implies that the Bond 
number. Bo = pgR^/a, is small. Here g is the accelera- 
tion due to gravity. For instance, even under terrestrial 
conditions, an air bubble of radius R = 0.1 cm suspended 
in water is characterized by Bo « 0.1. Thus, the impact 
of gravity on the bubble results in a relative surface dis- 
tortion on the order of only 10%. We are interested in 
smaller bubbles so that the distortion effects are insignif- 
icant. 

Because of symmetry, we use the spherical coordinates 
r, "d, a with the origin in the center of the bubble. As we 
have announced in Sec. IH we will demonstrate that the 
shape oscillations are able to interact with the breath- 
ing mode even within the linear approximation. Be- 
cause in this approximation the nonaxisymmetric modes 
of the shape oscillations do not interact with the breath- 
ing mode (see argumentation below in this section), we 
eventually restrict our consideration to the axisymmetric 
problem. 

We now formulate the dimensionless governing equa- 
tions and boundary conditions by measuring the dis- 
tance, time, velocity potential, deviations of pressure and 
the bubble surface in the scales of R, ^/pRF/a, t\J aRj p, 
ea/ R, and ei?, respectively. As one sees, the small pa- 
rameter e has the meaning of the ratio of the amplitude 
of the surface oscillation to the equilibrium radius R of 
the bubble. Small oscillations of the inviscid incompress- 
ible ambient are governed by the Bernoulli equation and 
the condition of incompressibility: 



(4) 



The solid surface = 7r/2 is impermeable for the liq- 
uid: 



dip 



0. 



(5) 



At the free surface, governed by the equation r = 1 + 
e t), we prescribe the kinematic and the dynamic con- 



dt 

^,3 



dip 
dr ' 
1 



P 
d 



Pg 



in 



2)C, 



(6) 
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where p is the velocity potential, p and Pg are the pul- 
sation parts of the pressure in the liquid and gas phases, 
respectively. These pulsation fields describe the pressure 
deviation from their equilibrium values: i~^PgQ, Pgo = 
PgR/a in the gas and {pgo — 2) in the liquid. 

Next we define the spatially uniform oscillations of the 
gas pressure Pg. Because the pulsations of the gas pres- 
sure in the bubble are spatially homogeneous and the 
dissipative processes are assumed to be negligible during 
a period of oscillation [recall the accepted restrictions 
([2])], we can apply the adiabatic law for perfect gas. As 
a result, we arrive at the condition 

[PgO+Pg{t)]V^t)^PgoVo^. 

Here 7 is the adiabatic exponent, Vq = |7r is the dimen- 
sionless volume of motionless bubble, and 



Vit) = Vo{l + 3e{0), (0 



1 

2^ 



Cds, 



where the angle brackets denote the space averaging over 
equilibrium surface S = 5(19) of the bubble. Thus, we 
obtain for the pressure pulsation in the gas phase: 



= -37p<,o (C) = -no (C) , 



(7) 



where we introduce parameter Hp = 3j PgR/a. Because 
Ho oc pgo, hereafter Hq is referred to as the dimensionless 
equilibrium pressure inside the bubble. Note that this 
parameter can be presented as a ratio of frequencies of 
the volume and shape oscillations squared, Hq — {(jJbTc)'^- 
Thus, one can clearly see that for finite values of Hq in- 
equalities H]) follow directly from requirements 0. 

We emphasize that for nonaxisymmetric modes (C) — 
0. As it follows from Eq. ([7|), the pressure inside the bub- 
ble for these modes does not change in time, Pg{t) = 0, so 
that the breathing mode cannot emerge. In other words, 
because in the linear problem the modes with different 
azimuthal numbers are noninteracting, no excitation of 
the volume oscillations is possible. Hence, to be able to 
focus on the announced interaction between the breath- 
ing mode and the shape oscillations hereafter we deal 
with the axisymmetric problem, {0 7^ 0. 

The formulation of the boundary value problem is com- 
pleted by prescribing the dynamics of the contact line, 
where we impose the Hocking condition ([3]) 



9C ^_,ac 

dt dd' 



(8) 



where the dimensionless number A = Ky/pR/a is the 
Hocking (also known as wetting) parameter. 
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The boundary value problem ([J])-® has been formu- 
lated in the linear approximation with respect to small 
e. Particularly, this approximation allows us to pose the 
boundary conditions ^ at the time-averaged, i.e., at the 
equilibrium, position of the surface, r — 1. The formu- 
lated governing equations and boundary conditions in- 
volves the two dimensionless parameters Xlg and A. The 
limiting case of high pressure, XIq ^ oo, refers to the 
consideration of incompressible gas. In this situation, the 
problem ([4])- ([8]) transforms to the one describing the nat- 
ural oscillations of an incompressible bubble immersed 
in a hquid^ (see a detailed discussion in Sec. IIVI) . Con- 
versely, for small values of Ho one expects that the bubble 
collapses. Note that hereafter we assume the pressure in 
the ambient fluid to be positive and are not interested in 
consideration of cavitation. The limiting situations with 
respect to parameter A correspond either to the fixed con- 
tact line (A 0, the contact angle can change) or to the 
fixed contact angle (A — > oo, the motion of the contact 
line is allowed). As it has been mentioned before, apart 
from these particular situations the Hocking condition 
leads to energy dissipation near the contact line. For this 
reason, the natural oscillations are generally damped. 

Let us represent the decaying at the infinity solution 
to the Laplace equation ^ , which satisfies the imperme- 
ability condition in the form 



E 

n=0 



AnP2n (0) 



„2n+l 



(9) 



Here we introduce the variable 9 — cos i) and the Legen- 
dre polynomials Pn{0). 

Substitution of ansatz ([9|) into the Bernoulli equation 
@ and kinematic condition [the first relation in Eq. ^\ 
leads to expressions 



p = Re 

c = 



^2n+l 



Re 



2 AnP2n{d) 

ji=0 
oc 

.n=0 



(10a) 
(10b) 



Next we determine C from the dynamic condition [the 
second relation in Eq. ©] 



C = Re 



^ OO 

:.no (C) - E 



tO^AnP2nie) 

(2n- l)(2n + 2) 



C9 



and equate it with Eq. (jlObp . Accounting for a relation 
{0 = —Aq, we obtain the coefficients introduced in ([9]): 



^0 



C 



An 



(4n + l)P2«(0)C 



(n > 0). 



(11) 

Here, = (2n— l)(2n+l)(2n-f 2) are the eigenfrequen- 
cies of the shape oscillations of a spherical bubble. These 
frequencies refer to the even modes; IIq = Hq — 2 is the 



frequency of the volume oscillations of a spherical bub- 
ble in liquid. Here, the term —2 presents the well-known 
correction to the frequency of breathing mode caused by 
surface tension. 

After the substitution of (fTO|) and (|lip into condition 
, we arrive at the dispersion relation defining the spec- 
trum of the eigenfrequencies of a bubble: 



E 



anP2n{0) 



1 



A (12) 



(4n + l)P2»(0) 
(2n- l)(2n-|-2)' 

where a„ are the coefficients in expansion of 6 in the 
series of the even Legendre polynomials at ^ G [Oi 1]- 

Note that except for some particular situations (e.g., 
small and high values of A) the eigenfrequencies defined 
by Eq. (|12p are complex. Separating the real and imagi- 
nary parts, one can demonstrate that the imaginary part 
of the eigenfrequency is non- negative for any Hp > 2, 
i.e., the oscillations decay. At Hq « 2 and a finite value 
of the wetting parameter, one of the eigenfrequencies is 
determined by the relation 



uj = a(Ho - 2), 



(13) 



which means that too small pressure inside the bubble 
results in the monotonic instability of the bubble with 
respect to collapse. Thus, relation (fT3|) defines the sta- 
bility threshold against the fast (adiabatic) compression. 
Clearly, at positive external pressure the bubble becomes 
unstable for higher values of the gas pressure: Ho < 67 
(or, in dimensional units, Pg < 2a /R). 

Consider now the dispersion relation in the limiting 
cases. For the fixed contact angle, A 00, the eigen- 
modes of oscillation of a hemispherical bubble coincide 
with the corresponding even modes of the spherical bub- 
ble. In this case, the oscillation frequencies equal 
(fc = 0, 1, . . . ). For high but finite values of the wetting 
parameter the eigenfrequencies obey a relation 



UJ 



(k) 



Ilk 

A 



Ik^k 

~2}^ 



( 

oc 

E 02 

n = 0, 



\ 



lk_ 



J 

(14) 

defined for fc = 0, 1, . . . . Here, 7^ = -O|afcP2fc(0)/2 
(fc > 0), 7o = 1/2. As it can be seen, the oscillations 
are weakly damped, a correction to the eigenfrequency 
is proportional to A~^. Formula is identical to that 
for the incompressible bubble for positive fc, if the first 
term in the sum is vanishing, which corresponds to the 
limit — > 00. We indicate that result holds even for 
Ho < 2. In this case all the terms for fc = are imaginary 
and relation (jl4[) describes two branches in the spectrum 
(the stable and unstable): w'"' = ±i^/2 - Hq + 0(A"i). 

The frequency of oscillation of the bubble with the 
fixed contact line Wp is determined from the following 
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FIG. 2: Eigenfrequencies uj of the bubble oscillations as func- 
tions of pressure Ho for modes — 2. The real (a) and imagi- 
nary (b) parts of u axe plotted for A = 1 (full line) and A = 10 
(dashed line). 
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FIG. 3: Eigenfrequencies u) of the bubble oscillations as func- 
tions of wetting parameter A for modes — 3. The real (a) 
and imaginary (b) parts of uj are plotted for Do = 30. 



real expression: 



IT. 

n=l 



.P2„(0) 



02 



= 0. 



(15) 



In the case of A <C 1 (slight shp of the contact hne) the 
results are quahtatively similar to those of incompressible 
liquids^i^l the decay rate is proportional to the small 
wetting parameter and a correction to the frequency ujp 
is proportional to A^. 

However, despite the qualitative similarity with the os- 
cillations of the hemispherical drop studied in Refs. [26ll27l . 
there arises a number of peculiar effects for the compress- 
ible bubble. The appearance of an additional branch in 
the spectrum of natural oscillation, governed by the pa- 
rameter XIq, leads to a nontrivial effect: the volume and 
the shape oscillations start to interact. Thus, if for a 
spherical bubble of the same radius the frequency Oq of 
the volume oscillations is close to any of the frequencies 
rik of the shape oscillations (for an even mode) 

then one of the eigenfrequencies wofc for the hemispherical 
bubble is in between these frequencies: 



70 Oa 



7/cOn 



70 + lk 70 + 7fc 



(16) 



A complex correction to the frequency defined by ([TC 
is proportional to S^. For this reason, the decay is weak 
for this mode irrespective of the wetting parameter and 
in the "resonant" situation {S — 0) is absent at all. With 
the accuracy up to the terms of order only two co- 
efficients, Aq and Ak, are nonvanishing in sums ^ and 
(fTU)) . Hence, this mode is a superposition of two kinds 
of motion: the radial pulsations and the fc-th mode of 
the shape oscillations of the bubble with the fixed con- 
tact angle. These oscillations occur in the antiphase and 
their relative amplitudes are such that the contact line 
remains motionless. 

For arbitrary values of the governing parameters, 
Eq. (fTi|) was solved numerically. Using the secant 



method we retained up to 200 terms. Careful analysis 
has shown that for most situations retention of only 20 
terms is enough to ensure convergence. 

In Fig. [2] we present the frequency of natural oscilla- 
tions as a function of the pressure in the bubble. At 
A = 10 the real part of the frequencies of the first modes 
almost coincide with the values Qqi ^i, O2. The imagi- 
nary part of the frequency is small and is described by for- 
mula (|14p . A distinction exists only in the case Oq ~ O^ : 
the imaginary part of the frequency drastically decreases 
for the volume oscillations and increases for the shape os- 
cillations. At moderate values of the wetting parameter 
(A = 1), with the growth of Hq the spectrum of the eigen- 
frequencies rearranges. In this case the damping time of 
oscillations is comparable with their period; the imagi- 
nary part of the frequency as a function of Hq is given in 
Fig. El^b). As has been indicated before, near the "res- 
onance" Oo ~ Ofc the imaginary part of the frequency 
for one of the modes tends to zero according to the law 
cj, ~ (5^ cx (Oo - Ofc)2. 

The dependence of eigenfrequencies on the Hocking 
parameter at a fixed gas pressure in the bubble can be 
found in Fig. [S] For modes 2 and 3, the behavior resem- 
bles the case of the incompressible liquidi^i^l the real 
part decreases with the growth of A, the imaginary part 
has a maximum at a finite value of the wetting parame- 
ter, turning to zero at small and high A. An interesting 
behavior of the eigenfrequencies can be observed for the 
two lowest modes, when the complex frequencies become 
close, see Fig. ID Note that an insignificant change of Hq 




FIG. 4: Rearrangement of the two lowest modes and 1 un- 
der small variation of pressure Hq. Eigenfrequencies ujr as 
functions of wetting parameter A, plotted for Ho = 30 and 
Ho = 25, cf. Fig. Ha). 
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is able to qualitatively rearrange the dependence uJr{X). 
We indicate that the exact coincidence of complex decay 
rates of the two modes resulting in such a rearrangement 
of the spectrum is possible only in a discrete number of 
points (A, Ho), which obey the complex- valued equation 
cc;W(A,no) =cc;(")(A,no). 

As the final point of this Section, we point out to what 
extent the presence of the substrate influences daniping 
of a bubble oscillation. Generally (see, e.g., Refs. 
oscillation damping of a spherical bubble may be caused 
by viscosity, thermal diffusion, and radiation of acoustic 
wave. Although all these phenomena are relevant for the 
hemispherical bubble, they can be remarkably changed 
by the substrate. 

As we stressed in Sec. HI understanding the impact 
of viscosity on the moving contact line is a notoriously 
complicated problem. However, an important estimation 
can be drawn from the problem with the contact line 
pinned,— A = 0. Accor ding t o this analysis, the dimen- 
sionless decay rate is \/ Tc/tv against conventional Tc/tv 
for a spherical bubble away from the solid surface. 

A similar situation is expected for the damping caused 
by heat diffusion. Note that this mechanism may be im- 
portant for both gas and fluid. Moreover, the ratios of 
thermal conductivities of the three media (gas, liquid, 
and solid) become important. If the conductivity of the 
solid is not negligibly small, the thermal boundary layer 
near the solid surface is developed. The decay rate is 
given by y/cubTt, whereas for a bubble in the absence of 
the solid surface one ends up with the usual value of UbTt . 

Acoustic irradiation has nondissipative origin and 
therefore the presence of the substrate is insignificant. 
Hence, the conventional condition ujbTa ^ 1 as in ([2]) is 
sufficient to neglect this phenomenon. 

As we see, the presence of the solid surface leads to the 
development of the boundary layers, which results in the 
faster oscillation damping compared to the case of a bub- 
ble away from the substrate. These changes, however, do 
not change restrictions ([2]). Thus, the criteria that allow 
us to neglect viscous and heat dissipation and acoustic 
irradiation loss remain conventional. 



III. FORCED OSCILLATIONS 

Consider the behavior of a gas bubble sitting on a plane 
solid substrate, as sketched in Fig. [S] We neglect gravity 
as before and now look at the problem of forced oscilla- 
tions. Assume that the substrate performs transversal vi- 
brations (with respect to its plane) with the amplitude a 
and the frequency uj; the substrate velocity in an inertial 
reference frame is awsinwt. We consider the frequency 
of vibrations to be comparable with both the frequen- 
cies of the shape and volume oscillations. We stress that 
spatially uniform pulsations are of crucial importance for 
the consideration of the bubble. This is in contrast to the 
problems addressing an incompressible drop,^^'^^ where 




FIG. 5: The geometry of the forced oscillations problem. 



the system is dominantly governed by the pressure differ- 
ence in the liquid. In our system, to define the uniform 
part of the pressure field one has to specify an additional 
condition away from the bubble. From the experimental 
point of view, the most convenient way to overcome this 
difficulty is to attach the bubble to the bottom of a liquid 
layer of depth H (large compared to the size of bubble) 
with the free surface," see Fig.[S] For a spherical bubble, 
a similar way of treating this peculiarity of the pressure 
field has been applied before, ^■'^^ where the bubble has 
been immersed in a vibrated column of liquid with a free 
surface and the impact of vibrations on the spherically 
symmetric oscillation mode has been studied. 

In the reference frame moving together with the sub- 
strate the small oscillations of the bubble are governed 
by the following dimensionless equations and boundary 
conditions: 



P 

r 

z 



_J_d(p 

TT dip 



1 



= 0, 



2 

dt dr 
h: p^O, 

TT dC 
2 ■ 'dt 



cos fit, V> = 0, (17a) 
(17b) 



1, ^ 



^^ = ^, 172p + Ho(C) = (V^ + 2)C(17c) 

(17d) 
(17e) 



A 



'di9' 



The problem (|17p has been nondimensionalized using 
the scales aH/R, paco'^H, aH a / pR^ for the deviation 
of the bubble surface from its equilibrium form, the pres- 
sure, and the velocity potential, respectively; for the dis- 
tance and time we use the same scales as before, in prob- 
lem (|4])-([8]). We note that the oscillations of the bubble 
surface can be considered to be small only provided that 



a < 



H 



(18) 



As we see, the restriction imposed on the amplitude a is 
much stricter than in the situation of an incompressible 
drop [a <^ R). This restriction, however, can be weak- 
ened for the weakly compressible bubble, see Sec. HVl 

Among the dimensionless parameters Hp and A, de- 
fined earlier, there appear two more parameters in prob- 
lem (|17p : the dimensionless frequency of the sub- 
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strate oscillation, which is related to the Weber num- 
ber Q'^ ~ puj^R^/a, and the relative width of the layer 
h = H/R. As we have assumed above, the last parameter 
is high: 



/i> 1. 



(19) 



This inequality allows us to neglect the term proportional 
to z/h in the Bernoulli equation (jl7ap . and apply an 
approximation. We replace the exact condition (jl7dp for 
the pressure at the free surface with a requirement 



OO, (/3 ^ 0, 



(20) 



which ensures the disturbance decay far from the bubble. 

Indeed, in the absence of the bubble the vibrations of 
the layer cause pulsations of the pressure in the liquid: 
Po = (1 — z/h) cosLut. We note that near the bottom, 
where z h, the inertial force is negligible, therefore 
these pulsations can be considered as spatially uniform. 
Thus, far away from the bubble but close to the solid 
surface, the time oscillations of the pressure are spatially 
uniform, which causes the volume (and hence also the 
shape) oscillations of the bubble. Surely, the wave scat- 
tered by the bubble does not satisfy condition (jl7dp pre- 
cisely. This solution, however, can be easily corrected 
by the method of images. Such a procedure introduces 
an error of order in the condition for the pressure 
at the bubble surface (|17c[) . which is of the same order 
compared with the neglected inertial term. 

The solution to the Laplace equation for the velocity 
potential that decays at infinity and satisfies the imper- 
meability condition (jl7b[) along with the fields of pressure 
and surface deviation, can be represented as follows: 




Tl = 



iSlt 



n=0 



(21a) 
(21b) 
(21c) 



Applying the dynamic boundary condition, we obtain the 
expansion coefficients 

n^ + C (4n+l)P2n(0)g , 

= 7^9 = 7^ 7^ (n > 0). 



(22) 

Here, the complex constant C is found from the Hocking 
condition (I17el): 



n=l 



a„P2»(0) 




As it can be seen, arg(yl„) is nonzero and identical for any 
n > 0, but different from arg(ylo). This indicates that the 



bubble oscillations is a superposition of the standing wave 
(shape oscillations) and the volume oscillations, which 
have phase shifts relative to each other and relative to 
the substrate vibrations. 

Further, one can show^^'^^ that the frequencies flk are 
not resonant at any value of A. At = flk the mo- 
tion of the bubble represents a combination of the radial 
(volume) oscillation and the k-th mode of the shape os- 
cillations (both these oscillations are in phase with the 
vibration motion of the substrate): 



c 



(P2k{0) 



-l]cosnkt {k>0), (24) 



which means that the contact line remains motionless. 

At a frequency close to fla, we put — flo{l + Sq), 
where the frequency mismatch Sq ^ 1, and obtain from 



a«P2„(0) 



$12 
"0 




(47i+l)P2„(0)r!§ 



$12 



{n > 0), 



C 



1 - Sonf^ % 



2 \ ^ 



a„P2«(0) 



n=l " 




Performing evaluation for — Q{6 — 0) and accounting 
for a relation X^J^o ^nP2n{^) — 0, we end up with 



Co ~ Re [ai^oe^^^"*] = -Ar^osinr^ot- 



(25) 



This result clearly indicates that at a finite value of A 
and a vibration frequency close to the frequency of the 
breathing mode, the bubble oscillates with a finite am- 
plitude. 

Next, in the limiting case of the fixed contact angle 
(A 3> 1) we arrive at an obvious conclusion: at any fre- 
quency the bubble performs radial oscillations with an 
amplitude 



A. = 



(26) 



Other coefficients in (|2T|) are small (of order A~^) because 
the contact line only weakly interacts with the substrate. 
If, however, the frequency $1 of the external force is close 
to the frequency Jl^ of the k-th mode of the shape os- 
cillations, then the resonant amplification of this mode 
occurs: 



c 



cos $lfet - 

7fc 



P2fc(0) 



(27) 



, tan/3fc = 



{Uk - 
Ik 



where 5o = — fi^). Thus, close to the resonant 

frequency (for the fixed contact angle) the bubble motion 
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FIG. 6: Amplitude-frequency response for the contact line 
^1 (a) and volume (b) oscillations evaluated for A = 1 and 
pressure values IIo = 50 (full line) and Ho = 20 (dashed line) . 



consists of a superposition of the radial oscillation and a 
standing wave (one mode in the expansion), which have 
a relative phase shift. Exactly at the point of resonance 
the solution is given by formula (|24p. 

For large A, in the situation where the frequency of 
oscillations is close to the eigenfrequency of the volume 
oscillations, f2 « Oqj a resonant amplification of the ra- 
dial pulsations takes place: 



c 



ao 



ao cos(rio< + /3o), 



(28) 



tan/3o — 



70 



Note that exactly at the point of resonance, = Jig, the 
obtained result is reduced to Eq. pS)) . which means 
that the two asymptotics match in the overlapping range 
of parameters: 57 « JIq and A 1. On the other hand, 
the amplitudes in ((28|) and relation ([26| become the same 
for Slo > 1^7 - > 

We also indicate that resonant solution (|28|) is char- 
acterized by oscillations with the amplitude oq = 0(A), 
which is much higher than the resonant amplitude Ofe at 
a frequency of the shape oscillations. As it follows from 
the resonant solution ((27)) . the amplitude ~ 0(1) and 
hence ao/ak ~ A ^ 1. 

We now consider the bubble dynamics when the three 
frequencies are close: « fig ~ ^fc • For arbitrary values 
of A we obtain 



C = A 1- 



P2k{0) 
P2k{0) 



COS ujkt, 



(29) 



i.e., the bubble oscillations are in phase with the sub- 
strate vibrations and the contact line is motionless. The 
amplitude of oscillations read 



A = 



7o7fc 



(30) 



where wofc is the eigenfrequency of oscillations for the 
regime flo ~ fife, described by Eq. (|16p . A more rigorous 
analysis indicates that the phase of the resonant oscil- 
lations is shifted by 7r/2 with respect to the substrate 
motion, and their amplitude is 



^k / 70 + Ik 



70 



^2 



(31) 
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FIG. 7: Absolute values of coefficients Ao (a) and Ai (b) as 
functions of near the double resonance, A = 1. Full and 
dashed lines correspond to Eqs. (I29|l and (|30|l at Ho — 15 
(5=1/12) and Ho = 20 (5=1/2), respectively. Squares (Ho = 
15) and circles (Ho = 20) present numerically obtained results 
according to Eqs. ([22|l and ((23|) . 



The resonant amplitude remains bounded even at the 
eigenfrequency wofc , except for the case of precise coinci- 
dence of all the three frequencies. Note, the fact that the 
amplitude of the bubble oscillations is inversely propor- 
tional to the frequency mismatch squared is in agree- 
ment with the result obtained in Sec. [TTl at ~ flk 
energy dissipation is proportional to S^. 

At arbitrary values of the governing parameters series 
PT|) were evaluated numerically. In Fig. [S] we present the 
amplitudes of the contact line oscillations Co = C{(^ = 0) 
and the radial pulsation as functions of frequency il for 
A = 1. Because of the dissipative processes at the con- 
tact line, the amplitude of resonant oscillations remains 
bounded. However, the interaction of the volume and the 
shape oscillations leads to considerable increase of the 
amplitude near the frequency defined by relation (|16p . 

As it follows from solution , the contact line is fixed 
at the frequencies Q coinciding with Qk, the eigenfre- 
quencies of the bubble oscillations with the fixed contact 
angle. At Q — ujp (recall that LUp is the eigenfrequency 
of a bubble with the fixed contact line) the amplitude of 
the contact line motion does not depend on the Hocking 
parameter: 



Cq{Q =^ LUp) 



(32) 



Note that for sufficiently small A, the function Co(fi) pos- 
sesses a local maximum at il — Up. Moreover, the local 
maximum is rather close to fl — LUp even at A = 1, see 
Fig. [6l For instance, the second maximum for Hq — 20 
(the dashed line) takes place at ^Imax = 11.89, whereas 



ujp = 11.94. 



As it becomes clear from Fig. [71 formulas ([^ and ([50]) 
work remarkably well at significant deviations from the 
"double resonance," e.g., the amplitude of the radial pul- 
sations at Ho = 20 ((5 = —1/2) is still in good agreement 
with (1201). 
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IV. OSCILLATIONS OF A WEAKLY 
COMPRESSIBLE BUBBLE 

We now turn to the consideration of the weakly com- 
pressible bubble, which implies high gas pressure, Hq « 
fig 3> 1. Compared to the situation analyzed in Sec. Illli 
the high pressure in the bubble means that the radial pul- 
sations as well as the induced shape oscillations become 
small. Mathematically, these facts are reflected by the 
smallness of C in (|21c[) : as it can be seen from Eq. ((22)) . 
all the coefficients A„ — > as fig ^ c». What is impor- 
tant, the condition Hq ^ 1 does not guarantee against 
the insignificance of the bubble compressibility. Indeed, 
the amplitude of the uniform part pg of pressure oscilla- 
tions is based on the width of the layer: pauj^H. This 
pressure contribution causes surface deviations propor- 
tional to IIq ^, it is of crucial importance for the bubble 
as a compressible object. Another source of bubble oscil- 
lations comes from the inertial force and induces bubble 
distortion independent of Ilg. The corresponding contri- 
bution in the pressure, pi„, is important for any bubble, 
does not matter compressible or not. This nonuniform 
part of pressure pin ~ pau^R, produced by the inertial 
force, is much smaller than the uniform part, because 
Pin/po — h^^ ^ 1. Thus, the dynamics of the weakly 
compressible bubble is determined by the competition of 
the two different factors: the weak compressibility itself 
(Ilg"^ <C 1) and the smallness of the inertial part of pres- 
sure (h^^ < 1). 

Formally, there appear two small parameters in the 
problem: h^^ and ■ Their ratio, the parameter that 
describes the impact of compressibility on the bubble os- 
cillations, is assumed to be finite. Note that restriction 
(dHI imposed on the amplitude in Sec. IIIII becomes much 
milder (see also Ref. 



h 



a < i?, a < ^R. (33) 

Thus, the solution to problem p7|) can be presented as 

(34a) 
(34b) 



(fi = iIq Vi + V2, 

p = cosilt + Ilg^pi + h^^p2, 



C = Hg^Cl 



(34c) 



Here, the fields ipi, pi, Ci describe the motion caused 
by the spatially uniform part of the pressure and ip2, P2, 

present the contribution induced by the inertial force. 
The latter contribution is related to the dynamics of an 
incompressible bubble. We also emphasize that the ac- 
cepted in Sec. [Ill condition Bo ^ 1 is sufficient to neglect 
gravity here, for the case of weakly compressible bubble. 
The smallness of the Bond number ensures that the cap- 
illary effects produced by the inertial force dominate over 
the contributions caused by gravity. 

The solution to the first problem, for the fields ipi, pi, 
Ci, is obtained from the results of Sec. IIIII in the limit 
Ilg ^ 1 (fig ^ 1). Writing down the fields of the veloc- 
ity potential ipi , pressure pi , and surface deviation Ci is 



terms of series ([211) , we obtain from ([22)1 the coefScients 
(4n+l)P2„(0)Ci 



{n > 0), (35a) 



^„-P2n(0) 1 , iX 



(35b) 



In Fig. [Sfa) we present the amplitude of the contact 
line oscillations Cio = Ci(^ = 0) as a function of fre- 
quency f2 for different A. It is clearly seen that at a fixed 
value of f2 this amplitude grows with the increase of A. 
This solution is closely related to the mode discussed in 
Sec. Uni Particularly, (ig = at = [see Eq. p4[) ]. 



At r2 = where is Up evaluated at the limit 

Ilg — > 00, the real amplitude of the contact line oscilla- 
tions is independent of A and equals il^ [ci. Eq. ([5^ at 
fig « Ilg > 1, recall the weight factor Ilg ^ in Eq. (|34cp ]. 
At A ^ 1 we conclude from Eqs. ([55)1 that Ci is small 
and Ci ~ —^0 ~ Note that this asymptotics works 

well even at A = 10 [see Fig. [HJa)], except for the close 
vicinity of "antiresonant" frequencies Q — flk, at which 
the contact line is motionless. 

The second problem, for (p2, P2, C2, is governed by the 
following equations and boundary conditions: 



P2 

d 

r 
r 



1 dip2 



2 -zcosVtt, VV2 = 0, (36a) 



2 

1 : 



1, -d 



= 0, 



dt 

TT _ dLp2 

8(2 _ dip2 
dt dr 
00 : (^2 = 0, 

TT ^ 5^ 

2 ■ dt 



(36b) 

f^'p2 - (V^ + 2) C2, (36c) 
(36d) 
(36e) 



dC2 



Problem ([36)1 describes the forced oscillations of an in- 
compressible bubble immersed in a liquid. Qualitatively, 
this problem is similar to that of the forced oscillations 
of a hemispherical drop,^^ Hence, the solution can be 
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FIG. 8: Amplitude-frequency response for the contact line of 
the weakly compressible bubble, (1 (a) and ^2 (b), evaluated 
for different values of A. The asymptotics of large A (left 
panel) is given by formula ^1 = —Q^. 
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represented as 

Lp2 = Re 

P2 = Re 



C2 



-Re 



^(2n+l)B„P2„(e)e' 



n=l 



(37a) 
(37b) 
(37c) 



where the expansion coefficients read 

. (4n + l)P2«(0)C2+a„ 



C2 = 



00 (2rt+l)a„P2„(0) 



Tl=l 



02 Q:„P2n(0) _ 1 1 iA 

" Z^n=i aj!-n^ 2 n 



(n > 0), (38a) 
(38b) 



The dependence of C20 = C2(^' = 0) on fl for solution 
(I57)) . ([55)1 is plotted in Fig. [HKb). As before, we observe 
a significant increase of the amplitude of the contact line 
oscillations with the growth of A. Note that at its own 
antiresonant frequencies, which are the zeros of the nu- 
merator of C2, we have (20 = irrespective of A. In 
particular, every standing wave, which is referred by in- 
dex n, has its own phase shift, i.e., solution p7|) with ((38|) 
describes traveling waves propagating along the bubble 
surface. Although not exactly the same, these phenom- 
ena related to the antiresonant frequencies and traveling 
waves are rather similar to those found in Ref. [26l . 

Having discussed separate contributions in Eqs. ([M)) . 
we can make the final note about the full solution. The 
expansion coefficients of this solution in terms of series in 
spherical harmonics are a superposition h~^Bn + IlQ^An, 
where both the contributions are small. The second term 
becomes negligible compared with the first one provided 
that 

pLJ^RH<^Pg. 

This requirement ensures that the compressibility of a 
bubble is negligible and incompressible approximation is 
valid. We note that a similar inequality has been applied 
in Ref. 32 in the context of bubbly media. 



V. CONCLUSIONS 

We have investigated natural and forced oscillations 
of a compressible hemispherical bubble put upon a solid 
substrate. The contact line motion has been taken into 
account by applying the Hocking boundary condition. 
We have proven that the linear shape and volume oscil- 
lations demonstrate interaction, which is the main qual- 
itative result of our study. Having performed detailed 
analysis, we have found out two important features ac- 
companying this interaction. First, we have detected a 



double resonance, where independent of the Hocking pa- 
rameter an unbounded growth of the amplitude occurs. 
Second, we have figured out the general condition that 
can be used to neglect the bubble compressibility. This 
requirement turns out to be stricter than it might be 
expected. Finally, although our analysis bases on the 
assumption of adiabatic oscillations, we indicate below 
how the obtained results can be used for an arbitrary 
polytropic process in the gas. 

We have focused on the natural oscillations and ana- 
lyzed the eigenfrequency spectrum. Particularly, we have 
addressed the question as to how wetting at the con- 
tact line, which is governed by the Hocking parameter 
A, influences the oscillation damping. We have neglected 
viscous and heat dissipation as well as acoustic irradia- 
tion loss. However, because the applied Hocking condi- 
tion includes its own energy dissipation mechanism, the 
eigenoscillations are generally damped. The dependence 
of eigenfrequencies and the decay rate on A are qualita- 
tively similar to the case of a drop^^ or an incompressible 
bubble Here, the decay rate is maximal for A = 0(1) 
and tends to zero at the limits of the fixed contact line 
(A — > 0) and the fixed contact angle (A 00). In the 
latter case, the eigefrequencies of the hemispherical bub- 
ble refer directly to those of the even eigenmodes for a 
spherical bubble of the same radius: the breathing mode 
frequency 0§ = Hq — 2 and the frequencies of the shape 
oscillations fll = {2k - l)(2fc -I- l)(2fc + 2), fc > 0. 

However, the compressibility of the bubble leads to a 
number of peculiar effects. Generally, the eigenfrequency 
spectrum becomes dependent on an additional parame- 
ter, dimensionless pressure in the bubble Hg. This de- 
pendence results in nontrivial interaction of the volume 
and shape oscillation and newly found rearrangement of 
branches in the spectrum, see Fig. U) As a particular 
consequence of the rearrangement, the eigenfrequencies 
of the compressible bubble are able to not only mono- 
tonically decrease with A, as in the case of the drop or 
incompressible bubble, but also monotonically grow. Of 
special attention is the "resonant" case when Oq ~ ^k, 
characterized by weak dissipation. At the exact coinci- 
dence of these frequencies the contact line is motionless 
and there is no dissipation irrespective of A. The bubble 
surface dynamics corresponds to a superposition of the 
fc-th mode of the shape oscillations and the antiphase 
radial pulsation. 

We have considered normal vibrations of the substrate 
and studied the problem of the forced oscillations. In this 
situation the main role is played by the spatially uniform 
pulsations of the pressure field, which cause the volume 
oscillations of the bubble. Through the interaction with 
the substrate via the moving contact line these volume 
oscillations induce the shape oscillations. The performed 
analysis of the forced oscillations has shown resonance 
phenomena to exist. Particularly, for weak dissipation we 
have obtained analytical expressions for the oscillation 
amplitudes valid close to the resonance. We have also 
found out the double resonance, f2 « fio ~ i^fc, where 
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is the frequency of substrate vibration. As it follows from 
pip , in this oscillation regime the resonant amplitude 
Ares oc (r^Q — f2|)^^. The divergence of Ares as ft 
flo = follows directly from the problem of natural 
oscillations, where the specific case fio = ^fc predicts no 
damping of oscillations. 

We indicate that although we have applied the adia- 
batic law leading to Eq. ([7]), our analysis holds for an 
arbitrary polytropic process. In this case, the adiabatic 
exponent 7, which enters Eq. ([7]) and the parameter IIq, 
should be now replaced with the polytropic exponent, 
m. This generalization makes our theory applicable to a 
wider range of bubble sizes and allows for the description 
of smaller bubbles, for which heat conductivity becomes 
the dominant dissipative effect.— »2, For instance, low fre- 
quency oscillations, uj <C XgR~^^ ^-re governed by the 
isothermal process, m = 1. We note, however, that this 
generalization is appropriate only for the forced oscilla- 
tions. For natural oscillations, any polytropic process is 
accompanied by intensive additional damping, which is 
different from the Hocking mechanism. Independent of 
A, this damping results in the complete attenuation of 
oscillations within a few periods. 

We have considered the special case of weakly com- 
pressible bubble and obtained the criterion identifying 
whether the bubble compressibility is insignificant. We 
have shown that the compressibility can be neglected 
only if the dimensionless pressure Ho in the bubble is 
large compared to the large /i, which is the ratio of 
the layer depth H to the averaged bubble radius R (see 



Fig. [5]). Our analysis allows us to draw a general conclu- 
sion about the impact of compressibility under the action 
of vibrations. It might be naively expected that com- 
pressibility effects become negligible at small frequencies 
LD in the sense that uj/lJc ^ 1, where Uc is a charac- 
teristic frequency of the volume oscillations. For bubble 
dynamics, uJc is the frequency of the breathing mode, 
tJb, for homogeneous fluid medicfc^^ it has the meaning of 
the acoustic frequency, ujc — c/H. However, the correct 
condition that does guarantee that the compressibility ef- 
fects can be neglected is significantly stricter and can be 
formulated as (uj/uJc)^ ^ /i- Here is a small parameter 
typical for a concrete physical situation. For instance, as 
we saw for the bubble dynamics ^ — h^^ <^ 1, whereas 
for thermoacoustic convection /i = fiQ <ti 1, where /3 is 
the thermal expansion coefficient and Q is the character- 
istic temperature difference. 
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